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SUMMARY 

Crackline-loaded edge- crack specimens are flat plate specimens which have a single 
crack notch extending normally from one edge and which are loaded in tension at posi- 
tions close to the intersection of the crack with that edge. These specimens are of inter- 
est in fracture -mechanics testing because they require comparatively little material and 
because there are some varieties for which the stress intensity factor K is almost inde- 
pendent of crack length over a considerable range. 

Stress intensity factors were determined for a variety of these specimens, all with 
straight boundaries, by boundary collocation of the Williams form of stress function. 

The boundary conditions were determined with the aid of another stress function due to 
Filon. The results presented are considered to be comprehensive enough for a wide vari- 
ety of applications of these specimens. 

Two complementary types of semifinite crackline-loaded specimen are treated in de- 
tail. The results for these specimens can be used graphically to obtain a useful prelim- 
inary estimate of the results for any finite specimen with straight boundaries. Boundary 
collocation analysis need then be conducted only on those cases which appear to be of 
most practical interest. 

The boundary collocation results were used to calculate values of the stress inten- 
1/2 

sity coefficient KBW ' /P for a crackline-loaded specimen of nonlinear contour. These 
values are in excellent agreement with the published experimental compliance data for 
this specimen configuration. 


INTRODUCTION 

The expression ” crackline-loaded edge- crack specimens” refers to members of 
that class of flat plate specimens which have a single crack notch extending normally 



from one edge and which are loaded by op- 
posed forces acting parallel to that edge, 
applied at positions close to the edge and 
also close to the crackline. The specimens 
discussed in this report have straight bound- 
aries with slopes Hp/e (= tan 2a), figure 1, 
varying from 0 to 0. 6. It should be noted 
that the convention used herein is that the 
dimensions a, e, and W, which are in the 
direction of the crackline, are measured 
from the line of action of the load P, not 
from the actual edge of the specimen. In a previous report (ref. 1) we presented stress 
intensity factors obtained by boundary collocation for one subclass of these specimens, 
namely, semifinite specimens (extending indefinitely beyond the crack tip) of zero slope. 
Detailed information concerning the use of crack-notch specimens for fracture- mechanics . 
tests of materials is contained in references 2 to 7. 

There are two aspects of crackline-loaded edge- crack specimens that are of partic- 
ular interest in connection with fracture-mechanics tests of materials and which led us to 
conduct a more comprehensive investigation of the relation between stress intensity fac- 
tors and crack length for such specimens. First, when used for plane strain crack 
toughness Kj c tests, specimens of this type can be designed so that they require less 
material per test than any type of remotely loaded specimen (ref. 2). In tests of tough 
materials of moderate yield strength, the required dimensions of specimens are a major 
limitation in determination of Kj c . The first such compact crackline-loaded specimen 
appears to have been that designed by M. J. Manjoine (ref. 3), which has since been ex- 
tensively used and studied (refs. 4 to 6). The current modifications of Manjoine’ s design 
are commonly called WOL (wedge- opening- loaded) specimens. 

The second attractive feature of crackline-loaded edge- crack specimens is that they 
can be shaped in such a way that the stress intensity per unit load is practically indepen- 
dent of crack length over a substantial range. This feature can be a valuable asset, for 
instance, in fatigue crack-propagation experiments or environmentally influenced crack- 
propagation experiments because the stress intensity can be controlled simply by con- 
trolling the load without regard to the crack length. The advantage of shaping specimens 
in this way was first brought to our attention by E. J. Ripling in 1964, and the results of 
his experimental studies with such specimens have since been reported (ref. 7). 

Although stress intensity calibrations of good accuracy have been published for a few 
special cases of crackline-loaded specimens (refs. 1, 4, and 7), we nevertheless consid- 
ered that a broad parametric study of the general specimen type would be valuable as a 
basis for judicious selection of specimen dimensions for specific purposes. The results 
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reported herein were obtained by the boundary col- 
location procedure that we have applied previously 
to other specimen configurations (refs. 1, 8, 9, 
and 10). This procedure is used to determine the 
coefficient of the first term of a Williams form of 
stress function (refs. 11 and 12), which is propor- 
tional to the stress intensity factor K. The neces- 
sary boundary conditions for the tapered specimens 
were determined by another stress function from 
Filon and Coker (ref. 13). 

It is possible to get a good general impression 
of the dependence of the stress intensity factor on 
the crack length for different values of the speci- 
men shape parameters (H , e, and W, or an equiv- 

__ alent set) by considering the two complementary 

^ types of semifinite specimen shown in figure 2. In 
one type, the distance W - a from the crack tip to 
the normal boundary is small compared with the 
crack length and with the distance to a transverse 
or oblique boundary. The crack-tip stress field is 
then dominated by the proximity of the normal 
boundary. It is convenient to refer to this case as the (W - a)-dominated type of semifi- 
nite specimen. There is an expression for a semifinite crack approaching the free edge 
of a half-plane given in the compilation by Paris and Sih (ref. 14, eq. (181)) which can be 
adapted to this case. The complementary type of semifinite specimen is that in which the 
normal boundary is far enough from the crack tip that its position has negligible effect on 
the crack-tip stress field. This type is referred to as the (W - a) -indifferent type of semi- 
finite specimen. This case was treated in reference 1 for specimens with zero taper (H 
constant), and it was shown that the dependence of KBH^^/P on a/H is linear. The 
present boundary collocation results show that the dependence of KBH^^/p on a/H^, 
where H a is the vertical height of the specimen arm at the crack tip, is also linear when 
the taper is not zero. Approximate treatment of these specimens as double cantilevers 
(refs. 7, 15, and 16) was not sufficiently accurate for our purpose, for reasons discussed 
in the appendix. 

Our purpose in considering the two types of semifinite specimens was to reduce to a 
minimum the number of boundary value computations that it was necessary to carry out. 

In conjunction with a moderate number of boundary collocation results to establish the 
necessary numerical coefficients, the semifinite relations can be used to sketch out the 
multidimensional field of the values of the stress intensity coefficient KBW^^/P with 
useful accuracy. Further computational effort could then be concentrated on those spec- 


lb) (W - a) - indifferent type. 

Figure Z - Complementary types of semifinite specimen. 
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imen shapes which are of most practical interest, as indicated by the preliminary results. 


SYMBOLS 


d~ , d, 


2n’ 2n-l 


e 

H 

H 

H 

H 

H 

K 

K 


w 


I 

Ic 


Units may be in any consistent system, since only dimensionless combinations are 
used in this report. 

a effective crack length 

B specimen thickness 

coefficients of Williams stress function 

distance from wedge tip to line of load application 

uniform depth of nontapered split arm 

depth of tapered split arm at crack tip a 

depth of tapered split arm at load line 

depth of tapered split arm at end boundary location 

depth of tapered split arm at location 6 

stress intensity factor of elastic stress field in the vicinity of the crack tip 

critical value of K at point of instability of crack extension in first or 
opening mode, a measure of plane strain crack toughness of material 

bending moment 

number of selected boundary stations for collocation solution 
load applied to both arms of specimen 
polar coordinates referred to crack tip 
specimen width 

cartesian coordinate system referred to crack tip 
half-angle of taper, 1/2 arctan (H^/e) 

distance measured from load line to new end boundary location 
normal direction to all boundaries 
cartesian coordinate system referred to wedge tip 
polar coordinate referred to crack tip 
polar coordinate referred to wedge tip 
distance from wedge apex to location C 


M 

m 

P 

r 

W 

x,y 

a 

6 

? 

v, 

e 

p 
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cr x stress component in x- direction 

Oy stress component in y-direction 

T^y shearing stress component 

cp polar coordinate referred to the wedge tip 

X Airy stress function 


BOUNDARY COLLOCATION PROCEDURE 


The method of analysis consists in finding a stress function that satisfies the bihar- 
monic equation v^y = 0 and also the boundary conditions at a finite number of stations 
along the boundary of a specimen, such as that shown in figure 1. The biharmonic equa- 
tion and the boundary conditions along the crack are satisfied by the Williams stress func- 
tion (refs. 11 and 12). Because of symmetry (fig. 1) the coefficients of the sine terms in 
the general stress function must be zero; hence, 


x( r >e) 
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The stresses in terms of y obtained by partial differentiation are as follows: 
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(a) Tapered-beam specimen boundary. 



(b) Force diagram about tapered-beam specimen, 
cantilever tip. 



(c) Rectangular-beam specimen boundary. 

Figure 3. - Geometric relations used in developing 
equations (3) and (4). 


The boundary collocation procedure consists in 
solving 2 m simultaneous algebraic equations corre- 
sponding to the known values of y and 9y/9^, where 
£ is the variable taken in the direction normal to the 
boundary (fig. 3(a)) at m selected stations along the 
boundary ABCD. Values for the first 2m coefficients 
of the Williams stress function are thus obtained when 
the remaining terms are neglected. Only the value of 
the first coefficient d^ is used for the present purpose 
since the stress intensity factor Kj is equal to -V 2 n dj 
as shown in reference 8 . 

For nonzero values of the angle 2a, the required 
values of y and its first derivative at the m selected 
boundary stations were obtained from the known solu- 
tion to the wedge problem with loads at the wedge tip 
(fig. 3(b)) as given in reference 13. The stress func- 
tion along boundary AB is formulated in the 77 , £ plane 
and is given in terms of the polar coordinates p and cp. 
This stress function may be transformed to rectangular 
coordinates x, y with origin at the crack tip by the 
following transformation (fig. 3(c)): 

x = p cos ( q ? + cp) - (e + a) = r cos 9 
y = p sin(ff + cp) = r sin 6 

where r and 6 are polar coordinates in the x, y 
plane. The equations for obtaining the boundary values 
are as follows: 
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Along AB, 




1. /- sin a sin <p + cos a cos <p \ + e /sin 2y - 2<p cos 2a _ A 

p P \2a + sin 2a 2a - sin 2a f 2p \2a cos a - sin 2a ) 


/ 3 2 ^ 

„ „ / -sin a . a cos o* - sin a cos a 

+ COS cp [ + 


1 2a + sin 2a 


2a - sin 2a 


/ 2 3 

(a sin a + cos a sin a . cos a 
- sin <p I — + 


2a + sin 2a 


2a - sin 2a . 


MU = Mil = MU caste + a) - 1 MU state * O) 

3? 3x dp p 3 cp 


Along BC, 


= (1 - 2 sin 2 a) - — 
P\P/ P 


9 W P > = -sin 2a 

3 ? 


Along CD, 


J-fX.) = (l - 2 sin 2 a) - -L 
P c \ p / P C 


Mx/P) = l(2^ = i(lZP) cos 2a - 1 d M Q sin 2a = 1 
3? 3x dp p dip 


V (3) 


where Pq is the distance from wedge apex to location C. 

The distance from the collocation boundary AB to the crack tip (the distance a - 6 
in fig. 3(a)) was varied for computational convenience. For 6 equal to zero, the resul- 
tant bending moment at boundary AB is zero, and the resulting shear force is equal to the 
applied load P. It is obvious that the boundary AB should not be taken too close to the 
crack tip on the basis of the St. Venant principle. 

A typical plot of the shear stress r and normal stress c r distributions along 

xy x* 

boundary AB (6 = 0) resulting from the stress function solution to the wedge problem is 
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Figured - Distribution of stresses o« and 
x^y at x= -a (fig. 3) based on stress runc- 
tidn solution to wedge problem where p - 1, 
fife - 0, e - 1. 5, B ■ 1, a - 4, W - 6, and 
Hp/e-0.20. 



Figure 5. - Distribution of stress and 
Tjjy at x = 6 - a (fig. 3) based on stress 
function solution to wedge problem where 
P -1, 6/e * 1. 07, e»L5, B - 1, a * 4, 

W = 6, and H p /e = 0.20. 


shown in figure 4. At this location, the resultant bending moment is equal to zero, and 
the shear load is equal to P. While the shear distribution at location AB is simi- 

lar to that resulting from simple beam theory, a small residual normal stress distribu- 
tion a occurs in the stress function solution. 

A typical plot of stress distribution when 6 in not equal to zero is shown in figure 5. 
The normal stress distribution cr x is very nearly the same as that obtained using simple 
beam theory, but the shear stress distribution is quite different. 

The stress intensity values obtained when 6 was taken to be different from zero 
were identical with those obtained for 6 equal to zero. 

For the special case of the rectangular specimen (Hp finite and 2 a equal to zero) 
as shown in figure 3(c), the equations for obtaining the boundary values are as follows: 
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Along AB, 


2L = 

P 


-6(a + x) jy J _ 


bh; 



-B-g^ =B ^Zg) = 
a? ax 




Along BC, 


Along CD, 


X - x + a 
P P 

p a <x/P) - F 9 (X/P ) „ Q 

a? ay 


( 4 ) 


2L = W 

P B 


B 9 (x/g) , p 9(X/P) . ! 
a? ax 


The boundary stations were equally allocated among the three boundaries; the spacing 
was uniform along each boundary. The collocation computations were carried out for 
each set of selected values of the primary variable a/Hp and the parameters 2 a and 
W/H . In general, the number of boundary stations m varied from 18 to 45. 


RESULTS AND DISCUSSION 
Preliminary Results for Semifinite Specimens 

As mentioned in the INTRODUCTION, the complete field of values of the stress inten- 
sity factor K can be derived with fair accuracy from a knowledge of the characteristics 
of the two semifinite specimens shown in figure 2. Equation (181) of reference 14 can be 
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kbh| /2 /p 


adapted to the (W - a)-dominated type of semifinite specimen by reformulating the relation 
in terms of a crack length a and a specimen width W. The result can then be expressed 
conveniently in terms of the dimensionless stress intensity coefficient KBW 1 / 2 /?: 


KBW 1 / 2 _ 0. 537 + 2. 17(1 + a/W)/(l - a/W) 
P ’ (1 - a / W)*/ 2 


There are several other forms of dimensionless stress intensity coefficient (for instance, 
KBHp/ 2 /p), but the selected form has the advantage that it is a single valued function of 
a/W for all values of the parameters Hp/e and W/e. Furthermore, it is appropriate 
for practical use since the measurement of W is more direct than the measurement of 
Hp or H a and has to be made in any case. 

To obtain quantitative expressions for the (W - a)-indifferent semifinite specimens of 
various slopes Hp/e, it was first necessary to obtain some appropriate boundary collo- 
cation results. The method proposed in references 7, 15, and 16 for deriving such ex- 
pressions by treating the specimen as a pair of cantilevers was not satisfactory (see ap- 


pendi ). When the necessary boundary collocation results are expressed in terms of 
KBH a ' 2 y/p and a/H a as in figure 6, it is seen that they are a very good fit to a family 
of straight lines with slopes A that depend upon Hp/e. These lines appear to have a 
common origin at a/H ^ -0. 7. We can offer no simple explanation of this apparently 

a 

simple finding; indeed, we suspect that the linear relations would not apply to actual 
specimens with sufficiently small a/H a . It must be remembered that the assumption of 

the St. Venant principle is implicit in 



-1 0 1 2 3 4 5 

a/H a 


Figure 6. - Boundary collocation results for (W - a) - indifferent 
semifinite specimens. Coefficients A of linear fitting expres- 
sions are also shown. 


the application of the boundary colloca- 
tion method, and this principle is vio- 
lated for small a/H„. Satisfactory re- 
suits could not be obtained by the bound- 
ary collocation procedure for (W - a)- 
indifferent semifinite specimens with 
Hp/e greater then 0. 4, presumably be- 
cause of the limitations of the computer 
program. It seemed reasonable, how- 
ever, to obtain approximate values of 
the coefficient A for Hp/e up to 0. 6 
by extrapolation of the values that were 
obtained. 

To show the relations for the two 
semifinite specimen types on a common 
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Figure 7. - Example of combined plot for two types of semifinite specimen. 


plot, the linear relations between KBH^y P and a/H a for the (W - a)-indifferent spec- 

imens are readily converted to relations between KBW— ^/P and a/W with parameters 
V e and either W/e or W/Hp. These latter relations are not linear except when 
Hp/e = 0. An example of a combined plot for the two types of semifinite specimen is 
shown in figure 7 (Hp/e = 0. 4). Individual boundary collocation results are also plotted 
on this figure to show how well they fit the curves. The agreement with the (W - a)- 
indifferent curves is good in the lower range of a/W, and there is fair agreement with 
the (W - a) -dominated curve when a/W = 0. 7. In the intermediate transition regions the 
boundary collocation results fall above both applicable curves. It is possible to construct 
a simple combination function which fits all the results fairly well, but there is little real 
point to doing so. 

By considering superimposed curves such as are shown in figure 7, we were able to 
select certain combinations of values of Hp/e and W/e which seemed to be of most po- 
tential interest. For each combination, we conducted boundary collocation analyses for 
values of a/W from 0.2 to 0. 7 at intervals of 0. 1. The results in terms of KBW^^/P 
are given in table I. 


Comparison With Experimental Results of Mostovoy et al. 

Experimental compliance measurement results are given in reference 7 for a speci- 
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TABLE I. - RESULTS OF BOUNDARY 
COLLOCATION COMPUTATIONS 


H/e W/e W/Hp 


W/H For a/W values of - 

Jr 

0.2 0.3 0.4 0.5 0.6 0.7 

KBW 1 / 2 /? 

1.25 5.31 6.86 9.29 13.44 21.52 

1. 5 5. 59 7. 10 9. 45 13 . 52 21. 54 

2 5.32 6.49 8.01 10.20 13.92 21.68 

3 7.65 9.35 11.34 13.53 16.64 23. 16 

4 10.20 13.00 15.76 18.57 22.35 26.81 

5 12.90 16.85 20.72 24.49 28.61 33.42 

>5 (a) (a) (a) (a) (a) (a) 


3 

30 

63.4 

70. 5 

74.0 

4 

40 

79.0 

85.2 

88.0 

5 

50 

92. 5 

98 

99 

2 

10 

17.65 

20.0 

21.3 

3 

15 

24.5 

26.7 

27.9 

3.2 

16 

27.0 

28.3 

29.0 

4 

20 

30.3 

32. 1 

32.8 

5 

25 

35.2 

36.2 

36.0 

3 

10 

14.7 

15.7 

16.1 

3.6 

12 

16.8 

17. 5 

17.7 

4.2 

14 

18.6 

19. 1 

19.2 

4. 8 

16 

20. 1 

20.4 

20.2 

5.4 

18 

21. 5 

21.6 

i 

21. 1 

6 

20 

22.9 

22.8 

21.9 

3 

7. 5 

10.2 

11.0 

11. 5 

4 

10 

12.4 

12.8 

12.9 

5 

12. 5 

14.0 

14.2 

14.2 

6 

15 

15.4 

15.5 

15.2 

8 

20 

18.2 

17.5 

16.8 

4 

8 


9.8 

10.3 

5 

10 


11.2 

11.0 

6 

12 


1Z 

11.6 

8 

16 


14 

12.8 




17. 8 19. 0 24. 0 

18.94 19.77 24.44 
19.85 20.46 24.82 


0.6 3 5 7 9.8 13.9 21.5 

4 6.67 7.7 10 14.0 21.6 

5 8.333 10.9 14.1 21.6 

6 10 10 11 14.2 21.7 

^se equation in fig. 6 or ref. 1. 



TABLE n. - CALCULATED VALUES OF KBW 1 / 2 /? 
FOR SPECIMEN OF REFERENCE 7 


a/W 

a / H a 

Hp/e 

A 

kbh^/ 2 /p 

KBW^/P 

Experimental result 
from ref. 7 

0.2 

1.04 

0. 53 

2.78 

4.85 

11.1 

10.9 

.3 

1.225 

.53 

2.78 

5.35 

10.8 

10.9 

.4 

1.36 

.45 

2.84 

5. 87 

10.8 

10.9 

. 5 

1. 50 

.45 

2.84 

6.26 

10.8 

10.9 


men which was contoured so that the stress intensity per unit load would be independent 
of crack length over a substantial range. It is of interest to compare these experimental 
results with the present results for (W - a)-indifferent specimens. The comparison is 
most conveniently made in terms of the stress intensity coefficient KBW^^/P. For the 
specimen of reference 7, the inferred value of KBW^/^/P is 10. 9 for all values of a/W 
between 0. 18 and 0. 55. The equation shown in figure 6 was used to calculate comparison 
values by the method of the previous section. To do this it was necessary to interpolate 
values of the coefficient A on the basis of values of H^/e calculated from the dimen- 
sions of the specimen of reference 7. (The exact dimensions are not given in the refer- 
ence; they were provided to us by the authors. ) The results are given in table n and 
agree very well with the experimental result. 


Final Results for Selected Specimens 

The precision of each result in table I is indicated in each case by the number of sig- 
nificant figures to which the value of KBW^^/P is given. The last digit is precise to 
within three units or less. Precision here refers to the variation of successive boundary 
collocation results with successive increase in the number of collocation stations beyond 
some minimum number. Some discretion has been excercised, however, in rejecting 
seemingly obviously disparate results which occurred occasionally in most runs. The 
precision is quite variable, and there are some blanks in the table where no useful re- 
sult could be obtained. The results for Hp/e = 0. 6 are the least satisfactory, and this 
slope appears to represent about the limit of capability of our computer program. 

For all values of Hp/e except zero, there are certain combinations of ranges of 
W/e and a/W in which KBW^^/P passes very gradually through a maximum or an 
inflexion. As mentioned in the INTRODUCTION, such ranges of almost constant 
KBW^^/P are of considerable interest for certain experimental purposes. Several ex- 
amples, selected from table I, are shown in figure 8. Clearly, there is a number of al- 
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ternative possibilities, and selection from 
among these for a particular purpose will de- 
pend on other considerations which are outside 
the scope of this report. One thing that should 
be realized, however, is that the proportions 
of all the ” constant- KBW^/1*" specimens in 
table I are such that the forward direction of 
crack propagation is unstable. Reference 7 
should be consulted for further information on 
this point and for guidance regarding the use 
of face grooves (or side notches) to circumvent 
this difficulty. 

For rectangular specimens (W/e = 0), 
KBW^/P increases more rapidly with a/W 
the higher the value of a/W, for all values of 
W/Hp. Again, practical considerations will 
determine the choice of proportions of a rec- 
tangular specimen for a given purpose. For 
K lc tests, there are two factors which 

strongly favor the choice of a relatively compact specimen (W/Hp about 2). These fac- 
tors are (a) economical use of material and (b) avoidance of the need for face grooving of 
the specimen to control the direction of propagation of the crack (ref. 2). 



Figure 8. - Selection of curves of KBW^/P and a/W 
for tapered specimens of such proportions that K is 
almost independent of a over a substantial range. 


CONCLUSION 

The boundary collocation results given in this report for specimens with straight 
boundaries appear to us to be sufficiently comprehensive for applications which we can 
forsee. Should more extensive results be needed, it is a simple matter to apply the re- 
sults for the two types of semifinite specimen to obtain a useful idea of the behavior of the 
appropriate dimensionless stress intensity coefficient in the region of interest. The 
boundary collocation method can then be applied to obtain more precise results for spe- 
cific cases. In this connection, the pattern of the precision of the boundary collocation 
results given herein provides an indication of where the method is most effective and 
where it is least effective. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 8, 1966, 

731-21-03-03-YBG0615. 
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APPENDIX - APPROXIMATE TREATMENT OF CRACKLINE-LOADED 
EDGE-CRACK SPECIMENS AS PAIRS OF CANTILEVER BEAMS 


Several treatments (refs. 7, 15, and 16) have been published in which the relative 
deflection of the arms of a (W - a)-indifferent type of semifinite specimen is calculated by 
the methods commonly described in texts on strength of materials. The result is used to 
determine the energy release rate with crack extension ^ from the Irwin-KLes relation 
(ref. 17): 


2 

= -E— — (2vB/P) (Al) 

2B 2 da 

where 2v is the relative displacement of the points of application of the load P. The 
stress intensity factor is given by 


K 2 = (A2) 

(for generalized plane stress, which is appropriate for a planar model). 

Since we were not able to obtain sufficiently accurate preliminary results by this 
method, it is useful to indicate its limitations. 

For small deflections, the local curvature of each arm of the rectangular specimen 
shown in figure 9 is 


d 2 y _ P(a - x) _ 12P(a - x) 
dx 2 EI EBH 3 


(A3) 


and the slope at any point is 



p 



Figure 9. - Rectangular (W - a) - Indifferent type of semifinite crackline -loaded edge- 
crack specimen. 
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ay = 

dx 



(a - x)dx + I Q (a, H) 


J *E-L-i.+ 

EBH 3 \ 2 



(A4) 


where I Q is the initial slope at 0 and is a function of a and of H, which cannot be de- 
termined by elementary methods because of the singular nature of the stress field at the 
crack tip. The deflection of the point of load application relative to the origin is, there- 
fore, 


y(x=a) = v = 


12P 

EBH 3 



and, hence, the compliance C of the specimen is 

* EBH’ 1 \' } / 


and 


dC 

da 



ebh 3 l 


+ 


_d_ 

da 




(A 5) 


(A6) 


(A7) 


therefore, from equations (Al) and (A2), 


K 2 B 2 


EB dC 
2 da 


12 

3 


H 


a 2 + — (I o a) 
da ° 


(A8) 


Since the term — (I a) is undetermined, this result is of little value unless a useful ap- 
da ° 

proximation can be found for the undetermined term. In this connection, it is useful to 
compare equation (A8) with the fitting equation to the boundary collocation results for the 
same type of specimen (ref. 1), as given in figure 6 of this report: 


KB H 1/2 ,3.46(^+0.t) 

P \H / 
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or 


K 2 B 2 _ U_ ^2 
P 2 H 3 


+ 1. 4 aH + 0. 5 H 2 ) 


(A 9) 


d 9 

so that -rjg (I Q a) is seen to be approximately equal to 1. 4 aH + 0. 5 H , which is by no 
means negligible compared with a 2 , except when a/H is large. This is not the case for 
specimens treated herein. 

Tapered specimens can be treated similarly by taking the angular bisector as the 
neutral axis of an arm. The analytical manipulation, however, is somewhat complicated. 
The outcome, as might be expected, is that the method is increasingly inaccurate with 
increasing specimen taper. 
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